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Abstract 

In this paper we use a simple model to explore the function of the gene Osteosarcoma-9. We 
are in particular interested in understanding the role of this gene as a potent anti-apoptotic factor. 
The theoretical description is constrained by experimental data from induction of apoptosis in cells 
where OS-9 is overexpressed. The data available suggest that OS-9 promotes cell viability and confers 
resistance to apoptosis, potentially implicating OS-9 in the survival of cancer cells. Three different 
apoptosis-inducing mechanisms were tested and are here modeled. A more complex and realistic model 
is also discussed. 
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I. INTRODUCTION 



Of the various Expressed Sequence Tags of unknown function which appear to be regulated 
when apoptosis is initiated one was identified as the gene Osteosarcoma-9 (OS-9 GI:580310), 
which is mapped to the chromosomal region 12ql3-ql5. This region is amplified in various 
human cancers, including sarcomas and osteosarcomas. The observation that none of the 
known genes in this region is amplified in all tested tumors suggests that other genes must exist 
that are commonly amplified in different cancer types. The functional analysis of the OS-9 gene 
is thus of interest for its potential anti-apoptotic role in the development of cancer. 

However, little functional data exist for OS-9. The identified sequence homology of OS- 
9 to both nucleolin and protein tyrosine phosphatase IB (PTP-IB) was found to entail no 
functional similarities 2|. Being driven from a TATA-less promoter [3|, however, suggests that 
OS-9 may be involved in cell-cycle regulated housekeeping or cell viability 2|. More recently, it 
was reported that OS-9 promotes oxygen-dependent degradation of hypoxia inducible factor la 
(HIF-la) via binding to both HIF-la and the HIF prolyl hydroxylases jj]. This is a surprising 
result since overexpression of HIF- la is characteristic of most human cancer and inhibition 
of HIF-1 impairs tumor growth. A role for OS-9 in the control of tumor progression as a 
tumor-suppressor gene has also been advanced 

In Ref. (3\ the first observational data set on a preliminary functional analysis of the product 
of the OS-9 gene was presented. The authors used overexpression of OS-9 isoform 2 cDNA 
in two murine FDC-Pl cell line populations, namely FDB-P1-17V (17V) and FDB-P1-17VIII 
(17VIII), and siRNA depletion of OS-9 specific mRNA in human umbilical vein endothelial cells 
(HUVEC) to study the role of OS-9 in induction of apoptosis by three mechanisms: interleukin 
3 (IL-3) deprivation of IL-3-dependent murine myeloid cells staurosporine treatment jsl, 
and tumor necrosis factor a (TNFa) treatment [iq|. The data on the number of viable and 
apoptotic cells in the two OS-9-overexpressed cell line populations were obtained from trypan 
blue exclusion analysis and from fiuoresence-activated cell sorter (FAGS) analysis after treat- 
ment with Annexin V and propidium iodide (PI). From the comparison with the cell population 
having the vector confering the Haemagglutinin (HA) gene only (negative control), the data 
suggest the following results: 1) promotion of cell viability upon induction of apoptosis by 
withdrawal of IL-3; 2) promotion of cell resistance to apoptosis upon induction of apoptosis by 
addition of staurosporine and TNFa; 3) concentration dependence of resistance to apoptosis by 
addition of staurosporine; 4) resistance to apoptosis, with HUVEC'S precipitating to apoptosis 
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upon blockage of the OS-9 function; 5) no effect on cell proliferation of any of the mechanisms. 

This paper proposes a mechanistic model for the OS-9 functional data presented in Ref. [6|], 
thus giving theoretical support to the interpretation of the experimental results by E. Vour- 
vouhaki et al. The model is general enough to avoid speculations about the particular dynamics 
associated with OS-9, but informative to the point of showing, in a quantitative manner, the 
impact of OS-9 in cell fate. In particular, such a model must reproduce induction of apopto- 
sis upon treatment with an apoptotic agent or upon withdrawal of a vital agent, and explain 
resistance to apoptosis and cell viability as a consequence of OS-9 overexpression. 

This paper is organized as follows. In section [TTl we introduce the minimal model for the 
kinetics of a population of murine FDC-Pl cells with OS-9 overexpressed when subject to 
three different mechanisms for induction of apoptosis, namely withdrawal of lL-3, addition 
of staurosporine and addition of TNFa. We analyse the biological system as a sequence of 
biochemical interactions which we then formulate into a mechanistic model encompassed by a 
system of kinetical equations. In section IIIII we optimize the parameters which describe the 
experimental results by fitting to the data the kinetical equations. For the optimal values for 
the parameters, the kinetical equations are expected to model the biological system within the 
observational errors and the limitations of a theoretical construct, as discussed in section HVl 



II. MODEL FOR OS-9 OVEREXPRESSION IN APOPTOSIS INDUCTION 

In this section we analyse the biological system described by the experimental data in Ref. 
as a sequence of biochemical interactions modeled by the overexpression of the OS-9 protein 
in two cell line populations of murine FDC-Pl cells, namely 17V and 17VIII. The level of 
OS-9 mRNA in 17VIII was confirmed by Real Time Polymerase Chain Reaction (PCR) to 
be approximately four times that of 17V [6]. Given the lack of knowledge on how the OS-9 
protein acts in the cell, the entailed biochemical interactions are reduced to the simplest sensible 
functional form. In the biological system, we identify two competing effects, namely induction 
of apoptosis on the one hand, and on the other hand promotion of cell viability and resistance 
to apoptosis due to OS-9 overexpression. 

Apoptosis is induced by the addition of an apoptotic agent, identified as a^, or by the 
withdrawal of a vital agent to the survival of the cells, identified as Va- Here we describe 
induction of apoptosis by a decay rate of viable cells in the presence of an apoptotic agent or 
absence of a vital agent. These rate parameters are characteristic of the apoptotic mechanism 
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and hence assumed to be dependent on the specific agent. 

The cell response to apoptosis induction will depend on the induction mechanism. Thus, 
addition of an apoptotic agent will provoke the cell to resist apoptosis, in which case we will 
be assessing the OS-9 action in resistance to apoptosis. Conversely, withdrawal of a vital agent 
will challenge the cell to strive for survival, in which case we will be assessing the action of 
OS-9 in promoting cell viability. 

These are the basic biochemical interactions which, together with the normal viable cell 
growth, must be contemplated by any model that attempts to describe the biological system 
being considered. We now proceed to formulate the simplest possible kinetical model containing 
just as many parameters as the data readily available allow us to constrain. An extension of such 
model for an arbitrary number of possible intermediary mechanical interceptions is presented 
as an illustration of the increased level of data required to constrain an increased level of 
interaction. 



A. Kinetical Model 



In our first attempt to model the 0S9 role, we collapse the biochemical interactions into 
a minimal set of kinetical equations. This first-order approximation avoids speculation and 
allows us to concentrate on the effects of OS-9 overexpression on promotion of viability and 
resistance to apoptosis. In the simplest plausible model, we have three variables: viable cells 
Ct,, apoptotic cells and dead cells c„. The evolution in time of the number of cells in each 
of these cell stages is dependent on the reaction of the cell population to the absence of vital 
agent Va or on the interaction with the apoptotic agent aa, and on the subsequent cell response 
derived from an OS-9 overexpression. 

The biochemical interactions can be analysed in fundamental phenomenological processes 
which we represent schematically by the following system of kinetical equations: 

Cy > 2ct,, (1) 
Cv ^ Cq,, (2) 

Ca ^ (3) 

Here Eqn. ([1]) corresponds to cell proliferation, Eqn. ([2]) to apoptosis induction by either addition 
of an apoptotic agent or withdrawal of a vital agent and Eqn. ([3]) to cell death. The reactions 
will now be described in detail. 
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Viable cells grow at a rate A. The absence of any effect of the OS-9 overexpression on cell 
proliferation, as concluded from the DNA content analysis, implies that A = constant. The 
effective cell growth, however, depends on the ratio of the number of viable cells to the carrying 
capacity c^^max)- As a first approximation, we assume that the cells proliferate effectively at a 
rate A(l — / c^^rnax)) ■ The choice of the linear functional form for the effective cell proliferation 
rate on the carrying capacity was motivated by the treatment of various aspects of cancer 
modeling discussed in References HQ. 

Equation ([2]) describes the transition from viable cells to apoptotic cells induced by either 
addition of an apoptotic agent or removal of a vital agent Va- For both apoptosis induction 
processes we use the same expression for the conversion rate: a{l — (3). The two parameters 
involved are used to segregate the effect of OS-9 overexpression. The parameter (3 is a, function of 
the OS-9 overexpresssion, (3 = /5([OS-9]), and takes values in the interval [0, 1]. The parameter 
a absorbs all other contributions to the conversion rate and depends on the specific agent 
leading to apoptosis. All experimental situations described in this paper use either one of these 
two induction processes; no experiment uses both. Although both apoptosis-inducing processes 
have the same rate expression, the interpretation given to the parameter j3 is different. In the 
case of vital agent withdrawal, the (3 parameter of Eqn. ([2]) measures promotion of cell viability 
modulated by the OS-9 overexpression. However, in the case of apoptotic agent addition, /3 
reflects resistance of a viable cell to becoming apoptotic. This process must thus take place 
upstream in the apoptotic signal pathway before the activation of the initiator caspases (caspase- 
8, -9 and -12) [13|. 

Cell death, as expressed in Eqn. ([3]), given by the expression 7(1 — 6). 

As done previously, this expression segregates the dependence on OS-9 overexpression. The 
parameter 5 is a function of the OS-9 overexpression, 6 = 5([OS-9]), and takes values in the 
interval [0, 1]. Apoptotic cells will enter the final stage of apoptosis at a rate 7. Correspondingly, 
the parameter 6 refiects resistance of an apoptotic cell to entering such stage characterized 
by the cell turning into an apoptotic body and eventually dying. This effect will take place 
downstream in the apoptotic signal pathway before the activation of the executioner caspases 
(caspase-3, -6 and -7) (iSj. 

It is important to emphasise the modulation produced by OS-9 overexpression in these 
reactions. This modulation is encapsulated in the parameters /5 = /3([OS-9]) and 6 = 5([OS-9]), 
which are such that (3,6 ^ [0, 1], where denotes the control (normal) expression and all other 
values reflect a contribution of OS-9 overexpression. 
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We now proceed to formulate the phenomenological model into a system of differential 
equations for the three variables. In the case of induced apoptosis by the addition of one 
apoptotic agent, the interaction can be broken down as follows: a) binding of to c^,; b) 
interception by the OS-9-codified protein of a signal pathway triggered from an effective binding; 
c) blocking of apoptosis induction in or blocking of necrosis induction in Cq, depending on 
the position in the signal pathway of the interception; d) survival or death of Ca- 

For a contact rate a, a molecules of the apoptotic agent will bind to a cell per unit time. 
If we assume that from each contact point a signal for apoptosis is always started, then ac„ 
will be the number of triggered apoptotic signals per unit time. Before starting the upstream 
apoptosis- inducing caspase chain reaction, (5[ac^ apoptotic signals interact with a molecule of 
the OS-9-codified protein per unit time, resulting in the interruption of the signal pathway. The 
remaining (1 — /5)(aCt,) signals will successfully induce apoptosis. In cells where apoptosis has 
already been triggered, and in the absence of OS-9 overexpression, the final stage of apoptosis is 
reached at a rate ^Ca- In the presence of OS-9 overexpression, however, for 5(7Ca) successful OS- 
9-led interruption of a downstream signal for necrosis, the final stage of apoptosis will instead 
be reached at 7Ca(l — 5). Combining the reactions described above, we find that the variation 
in time of the number of viable cells is given by 

^ = Ac, f 1 - - ac.(l - (4) 

Given the time scale of the experiment, cell decay is rendered negligible and thus unaccounted 
for. This is supported by the data recorded in the FACS diagrams in Ref . Q| , which suggest that 
necrosis is sourced by apoptosis rather than by decay. It follows that the number of apoptotic 
cells will vary in time as 

dca 

= ac„(l - /5) - 7Ca(l - 5), (5) 



and the number of dead cells as 

dc 



This is an example of a quasispecies-type system 



In the case of induced apoptosis by withdrawal of a vital agent, the interaction can be broken 
down in a similar way to the previous case, with the apoptosis signal now being driven through 
the mitochondria/cytochrome c pathway [7|. The induction of apoptosis can still be described 
by a decay mechanism of rate ac^. This means that a apoptoptic signals are triggered per unit 
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time in a cell subject to a vital agent-depleted medium. Of the total ac^ apoptotic signals 
generated, P{aCy) are intercepted by the OS-9 protein, which measures the rate of failure of 
apoptosis induction. 

Both mechanisms discussed here are independent and each comprises either addition of an 
apoptotic agent or withdrawal of a vital agent. Consequently, we would expect that each 
highlights an independent function of OS-9, respectively promotion of viability or resistance to 
apoptosis. 



B. Chemical Model 



The model described in the previous subsection can be regarded as a black box model which 
acts on a given initial state to reproduce the observed final state in a OS-9-dependent way. If 
the detailed modus operandus of the OS-9 action in the cell was known, then the system would 
be approximated by a multiple little black box model, with smaller black boxes describing inter- 
mediate reactions triggered by intermediary apoptotic agents and constrained by intermediary 
initial and final known states. For one intermediary apoptotic agent, a'^\ in the upstream 
pathway between the trigger of the exterior apoptotic agent on the membrane of the cell and 
the interception by the OS-9 protein which binds to the signal at the rate C(up{i), the variation 
in time of the number of viable cells will be given by 



dcy 
It 



Ac, 1 - 



C71 



\Cy\l- 



C71 



^v{max) _ 



— acy 



(2 - /3) - 1 



(7) 



for two intermediary upstream apoptotic agents, a^^^ and a^^\ with binding rates Q!up(i) and 
oi.up{2) respectively, we have that 



dCy 

'dt 



— acyaup{i)aup(^2) 



+ acy 



XCyll- 



^v{max) 

/3q;„p(i)Q;„p(2) + aup(i) (l — o:up{2)) + (l ~ 

Ctj 



— aCy 



0(up(l)<^up(2) (2 - /3) - 1 



^v(max) , 

and similarly for n intermediary chemical agents upstream in the apoptic pathway 

Cy 



dCy 



(2 - /3) n 



i) 



(8) 



(9) 



\ ^v{max) , 

To account for the OS-9-modulated interception of the downstream signal pathway we must 
allow for intermediary downstream necrotic agents, n'^\ binding to the OS-9 protein at rates 



adown{n)- It follows that 



dCg 

dt 



(2-/?)n««p( 



i) 



1 



i=l 



(2-5)n«d. 



'.own{i) 



1 



i=l 



(10) 



The role of each of these agents in the transmission of the apoptotic signal must be known and 
included in the dynamical system. To study the potential role of such a weakly known gene such 
as the OS-9, this model is therefore inadequate. As the logical extension of the kinetical model 
to a series of such smaller systems, this model could be applied to a multiple level interaction 
approach to biological systems with available multiple level data constraints jjjj] . It might also 
provide further insight into some unexpected results obtained with the kinetical model. We 
hope to pursue this in the future. 



III. EXPERIMENTAL AND THEORETICAL RESULTS 

In this section we proceed to use the system of first-order, coupled differential equations 
derived in the previous section to model the data reported in Ref. j^. This system is used 
to model the two different apoptotic inducing mechanisms, namely addition of an apoptotic 
agent or withdrawal of a vital agent. In this paper we shall only use what we called the kinetic 
model. We want to solve for the time evolution of the number of viable cells c„(t) and of 
apoptotic cells Caif) using the data in Ref. p] to constrain the parameters a, /5, 7 and 5 of the 
differential equations. In particular we are interested in the values for the parameters (3 and 5 
which measure the contribution of the OS-9 overexpression in the system. 

The data consist of the average of nine readings of the number of viable cells and the number 
of apoptotic cells obtained from trypan blue exclusion analysis and from FAGS analysis after 
treatment with Annexin V and PI, all experiments having been performed three times and 
each in triplicate. The FAGS analysis sorts cells that are 1) alive, 2) undergoing apoptosis, 
3) undergoing necrosis and that are 4) in between apoptosis and necrosis, thus allowing to 
discriminate alive, apoptotic and dead cells. The trypan blue exclusion analysis simply separates 
dead from alive cells, lacking the sophistication and sensitivity to discriminate viable cells from 
apoptotic cells. In a modeling study that considers all three of the cell stages, the data from 
the typan blue exclusion analysis are inappropriate for use and hence will not be included here. 

The solutions were determined numerically with mathematica (Wolfram Research). For 
the parameter estimation we implemented a sampling algorithm which calculates a cost function 
of a list of values for each parameter that we want to estimate. The cost function is defined as the 
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sum of the square of differences between the experimental points and the corresponding points 
of the solution to the kinetical system. The minimum of the cost function line, or surface, defines 
the parameter value which optimally fits the experimental data. The precision of the method 
is defined by the density of points in the search space. The optimal values thus determined 
for all the parameters were then used to plot the solutions for Cy{t) and Ca{t) which we present 
in this paper. We applied the algorithm to both the average values and the corresponding 
maxima and minima encompassed by the error bars. For each data set, we obtained from 
three to eight estimates for (and depending on) each parameter by fitting the equations to the 
average and to the extreme data points. The parameters obtained from the average data points 
were combined to produce the induced average solution (represented by a continuous line); the 
parameters obtained from the extreme data points were combined to produce the maximum 
and minimum solutions (represented by dashed lines siding the continuous line), which reflect 
the propagation of the measurement errors and thus define the induced error interval of all 
possible curves capable of optimally fitting the data within the experimental errors. All the 
curves show percentages normalized to the initial number of viable cells, c^,(t = 0) = 1 x 10^. 

From the experimental curves we can extrapolate the values for A and Cy(^rnax) adjusted to 
the system. Here we use A = 0.01 and Cv(^rnax) = 20 x 10^. From the data on the HA-only 
marked (negative control) population we estimate the parameters which do not depend on the 
OS-9 overexpression, namely a and 7. We proceed to apply them to the OS-9-overexpressed 
cell line populations to estimate the parameters which measure the overexpression of OS-9, 
namely (3 and 6. The data analysis of each apoptosis-inducing mechanism is introduced by the 
interpretation advanced in Ref. j^l based on a qualitative reading of the experimental results, 
which we then confront with the interpretation of the quantified, data-constrained parameters 
of the model. 

A. Withdrawal of IL-3 in FDC-Pl cells 

Here we analyse apoptosis induction by IL-3 deprivation of IL-3-dependent cells. The cor- 
responding apoptotic mechanism is that of withdrawal of a vital agent and is described by the 
system of equations (jlH6]), with /3 having the interpretation of OS-9 overexpression modulation 
in cell viability. The experimental results are collected in Figs. 2B and 2C of Ref. which we 
reproduce in Figs. [T]^a) and [2](a) respectively. 

First we solve for c„(t) given by Eqn. (jlj) and fit it to the number of viable cells shown in 
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Fig. [11(a). In the absence of OS-9 overexpression this equation has only one free parameter, 
namely a. Setting /5 = 0, we fit the solution to the curve for the number of viable cells in the 
HA-only marked cell population to find a = 0.057lo;oo}. In the presence of OS-9 overexpression 
we have additionally the parameter /5 to consider. We now use the results found for a to 
fit c^it) to the experimental curve for the number of viable cells in the OS-9-overexpressed 
cell hne populations and thus to find the corresponding /?. We obtain /^lyy = 0.17lo;o5 and 
Pirviii = 0.3llo;o4 for the 17V and the 17VIII cell line populations respectively. The resulting 
plots are depicted in Fig. [11(b). 

We then solve for Ca{t) given by Eqn. ([5|) and fit it to the number of apoptotic cells shown in 
Fig. [21(a) to determine the parameters 7 and 6. Using a = 0.057to;oQi and setting P,6 = 0, the 
fitting of the solution to the number of apoptotic cells in the negative control population gives 
7 = 0.30i[5:[jt Substituting now also for (3 the resuhs pnv = 0.17^^^ and pnviii = O.Slt^of^, 
the fitting of the solution to the number of apoptotic cells in the OS-9-overexpressed cell line 
populations gives 5nv = 0.099^o;o99 and Snviii = 0.086]'^o;o86 respectively. The resulting plots 
are depicted in Fig. [2]^b). 

In accordance with the experimental data, our mathematical model returns no significant 
difference between the number of apoptotic cells in the negative control population and that in 
the OS-9-overexpressed cell line populations, as illustrated by the shared pattern of the curves. 
We note that it captures rather closely the crossing point where the number of apoptotic cells 
in the OS-9-overexpressed cell populations is the same and the starting dominance of the 17V 
population over the 17VIII population is reversed. As anticipated by E. Vourvouhaki et al. 6|, 
this case illustrates promotion of cell viability (as shown by /? 7^ 0) and possibly downstream 
resistance to apoptosis (as shown hj 6 7^ 0). Our mathematical model allows to separate 
the effect of promotion of viability, as encapsulated in a non-vanishing j3, from the effect of 
downstream resistance to apoptosis, as encapsulated in a non-vanishing 6. 



B. Staurosporine treatment in FDC-Pl cells 

Here we analyse apoptosis induction by treatment with staurosporine. The apoptotic mech- 
anism is that of addition of an apoptotic agent and is described by the system of equations 
(|3H6]), with (3 having the interpretation of OS-9 overexpression modulation in apoptosis resis- 
tance. The experimental results are collected in Figs. 5B and 5C of Ref. jsl]. 

Two different concentrations of staurosporine were considered: l/iM and O.l/xM. The 
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data for [staurosporine]=l/iM are reproduced in Figs. [3](a) and ll](a); the data for 
[staurosporine]=0.1/iM are reproduced in Figs. [5](a) audita). We do not impose any addi- 
tional constrain in the system to account for the concentration change. In other words, we let 
the parameters a and /5 absorb the changes induced by the different concentrations. 



1. [Staurosporine]=liJ,M 

First we solve for Cy{t) for [staurosporine]=l/iM, given by Eqn. (jl]), and fit it to Fig. M^a) 
to find a = 0.20lo;o?^ from Cy in the negative control population, and Pnv = 0.18lo;o5^ and 
Pirviii = 0.24lo;o7 from the corresponding OS-9-overexpressed cell line populations. The results 
are plotted in Fig. [3](b). 

We then solve for Ca{t) for [staurosporine]=l/iM, as given by Eqn. ([5]), and fit it to Fig. Hl^a) 
in order to determine the parameters 7 and 5. Setting P,S = and using a = 0.20to;Q'j'^, we fit 
the solution to the experimental curve for the number of apoptotic cells in the negative control 
population and obtain 7 = 0.72lo;Jo- We proceed to fit the solution to the experimental number 
of apoptotic cells in the OS-9-overexpressed cell fine populations, using the values determined 
above for the other parameters, to find that 5nv = 0.19to;}g and Snviii = 0.27to;i3. The results 
are plotted in Fig. \Mjo)- 



2. [Staurosporine] =0.1 

First we solve for Cv{t), given by Eqn. (jlj), and fit the solution to the curves for depicted 
in Fig. [5]^a). Setting /3 = 0, we fit the solution to the experimental curve for in the negative 
control population and find a = 0.15l;o oo5. Fitting the solution to the experimental curves for 
Cy in the OS-9-overexpressed cell line populations with a = 0.15to;oo5 we find Pnv = 0.19to;Q5 
and Pirviii = 0.351°;°^. 

We then proceed to solve for Ca{t), given by Eqn. ([5]), and fit the solution to the experimental 
curves for Ca depicted in Fig. [6](a). Setting P,S = and a = OAbt^'^l^, we fit the solution to 
the experimental curve for Cq in the negative control population and find 7 = 0.20l!;o;o5^. This 
value is substantially lower than the one obtained for [staurosporine]=l/iM (7 = 0.721q;|o)- 
This difference emphasises the strong dependence of cell death rate with the apoptotic agent. 
Using now also the values for (3, we fit the solution to the experimental curve for Ca in the 17V 
and 17V111 cell fine populations to find both 6i7v = and Snviii = 0. These results are plotted 
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in Fig. [S](b) and Fig. [n](b). 

As anticipated by E. Vourvouhaki et al. j^, this case illustrates resistance to apoptosis, 
primarily upstream {(3 ^ 0) and possibly also downstream (5 7^ 0), dependent on the concen- 
tration of staurosporine. This is supported by the observed dependence of the apoptosis time 
scale on the concentration of the apoptotic agent. A higher concentration of staurosporine 
increases the probability of interaction with a viable cell, as encapsulated in the larger value of 

«(i/iM) = 0.20 > a(o.i/iAf) = 0.15. The probability of survival is further reduced by the ob- 
served decrease of upstream resistance to apoptosis with a higher concentration of staurosporine, 
which can be accounted for by the smaller value of /3, iirjyix^M) = 0.18 < /3i7y(o.i^A^) = 0.19 and 
(ivjviii{\ii.M) = 0.24 < /?i7v///(o.i/iM) = 0.35. The observed higher dying rate of apoptotic cells, 
on the other hand, can be accounted for by the larger decay rate 7(i^m) = 0.72 > 7(o.i/iM) = 0.20. 
The dependence of /3 and 7 with the concentration of the apoptotic agent was, however, not 
expected in the phenomenological framework of the kinetical model. 

C. TNFq treatment in FDC-Pl cells 

Here we analyse apoptosis induction by treatment with TNFa (40ng/ml=0.78yuM). Similarly 
to the previous case discussed, the apoptotic mechanism is that of addition of an apoptotic 
agent, which is described by the system of equations (j4H6]), with /5 having the interpretation of 
OS-9 overexpression modulation in apoptosis resistance. The experimental results are collected 
in Figs. 6B and 6C of Ref. which we reproduce in Figs. [Tl^a) and[8t^a) respectively. 

We want to solve for Ct,(t), given by Eqn. (jlj), and fit it to Fig. [Tl^a). Setting /3 = 0, we fit the 
solution to the experimental curve for c„ in the negative control population cell and find a = 
0.32+ 

aoi- Given a, we consider the case where P 0. Fitting the solution to the experimental 
curve for Cy in the OS-9-overexpressed cell line populations, we find Puv = 0.461q;o4 and 
f3^7Viii = 0.5Atom- The resuhing plots are shown in Fig. [7](b). 

We then proceed to solve for Ca{t), given by Eqn. ([5]), and to fit it to Fig. [8]^a). Using the 
values for a and /3 determined above, we fit the solution to the experimental curves for Ca in 
the negative control cell population and in the OS-9-overexpressed cell line populations, finding 
respectively 7 = 0.045to;Qjg and both 5i7v = and Suviii = 0. The resulting plots are shown 
in Fig.EKb). 

From the parameter estimation we conclude that this case illustrates upstream resistance to 
apoptosis {(3^0), in agreement with what was anticipated by E. Vourvouhaki et al. Q|. 
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IV. DISCUSSION 



In this paper we present a simple theoretical model to assess the possible effects of 
OS-9 overexpression upon induction of apoptosis. We formulated a mechanistic model for 
the reaction of OS-9-overexpressed cell populations to induction of apoptosis by different 
mechanisms and estimated the parameters of the model from the data set available on the 
possible implication of OS-9 in promotion of cell viability and resistance to apoptosis. We 
summarize the results in Table [B 
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TABLE I: Kinetical model parameters estimated for both cell line populations. 



From this Table we observe that the OS-9 overexpression- dependent parameters take non- 
vanishing values. Moreover, comparing the /3 values across all the apoptotic inducing mecha- 
nisms, we observe a strong dependence with the level of OS-9 overexpression in the two cell 
line populations 17V and 17VIII. These observations suggest that OS-9 has a possible role 
as an anti-apoptotic factor and that this role is dose-dependent, with the cell line population 
17VIII exhibiting a greater effect in both promoting cell viability and resisting to induction of 
apoptosis than 17V. 

A closer look at the S values reveals that they are very close to zero except for one situation, 
namely [staurosporine] = 1/iM for the 17VIII cell line population. We are persuaded to specu- 
late that the OS-9 overexpression does not have a significant effect on downstream modulation 
of apoptosis. The vanishing values found for the parameter 6 whenever presented without error 
bars are strictly vanishing, i.e. in the parameter interval [0, 1] the minimum of the cost function 
is at 5 = 0. In the cases where the induced average value for 6 is non-vanishing, the induced 
error bars can still contemplate a vanishing 6. A possible way of testing this effect would be to 
measure the OS-9 overexpression in the apoptotic phase only, as described by Eqn. ([3]). The 
number of apoptotic cells would be controlled, so that the conversion rate will only depend on 
7 and 6. 
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We observe that our model reflects the tendency observed in the data collected in Ref. 
and supports the interpretation therein advanced. There is quantitative evidence that OS-9 has 
a role in the regulation of the conversion rates for — > Cq and — c„, as illustrated by the 
non- vanishing values found for the OS-9-dependent parameters. Such evidence, however, is not 
determinant since the values are small and the errors are large (ranging from around 10-100% 
of the values estimated from the average data values). The results show that OS-9 contributes 
primarily to promotion of viability and upstream resistance to apoptosis, as reflected by /5 7^ 0, 
with little or no effect on downstream resistance, as reflected by the close to vanishing values for 
5. Furthermore, this function does not appear to be related to an increase in cell proliferation, 
as DNA content analysis shows that cells overexpressing 0S9 do not have a greater proportion 



of dividing cells than control populations [g]. 

Some results, however, remain to be fully explained. We collapsed in the parameter a all the 
intricacies of the intra-cellular, molecular processes involved in apoptosis induction and whose 
relation with OS-9 we had no data to constrain or infer. TNFa gives us the highest rate of 
apoptosis induction among the various other factors, followed by [staurosporine]= l/iM and 
[staurosporine]= O.l/xM. The withdrawal of IL-3, being a factor that promotes cell viability, can 
be seen indirectly to affect cell viability. Of the three apoptosis induction mechanisms described 
in the available data sets, only the staurosporine-mediated one allowed to discriminate the effect 
of the concentration of the apoptotic agent, [a^]. For that data set we observe a dependence of 
a with However, our model cannot predict the functional form of this dependence so as to 
reproduce the value A = 0.1251o;oq5 in the empirical expression 

O^jlf^M) ^ f — Y (11) 

relating the ten-fold increase of [oa] with a. Moreover, we observe a dependence of both (3 and 
7 with [aa] which is not a priori encompassed by the kinetical model, since these parameters 
measure the rate of cellular processes whose developement has, to our current understanding, 
no relation to the environmental conditions that triggered the preceding process. We also note 
that both /5 and 7 act towards a decrease of Cq, though in different ways. We venture that these 
observations might be evidence of the limits of the descriptive power of the kinetical model 
in face of the complexity of the biological system that it was conceived to describe as a first- 
order model. The binding rates aup{i) and adown{i) in the chemical model conceal intermediate 
phenomena at levels of interaction unaccounted for by the kinetical model. We believe that a 
higher-order model would be able to expose the sources of these observations, as well as allow 
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for possible dependencies with the apoptotic agent or the apoptosis inducing mechanism. 

Additional data from various experimental setups and from complementary measurements 
are, therefore, required to unambiguously determine the function of OS-9 in cell fate and in 
tumour growth, and test the model as a realistic theoretical description of the biological system. 
In particular, further experiments should attempt to identify the signalling pathway through 
which OS-9 exerts its function. Biochemical assays and western blot analysis will provide the 
means for the identification of the proteins that OS-9 binds to in order to transmit its potential 
signals. It is also essential that the function of OS-9 should be tested in in vivo models. The 
establishment of knockout mice is also important in order to assess the functional importance 
of OS-9 during development. Only then can we rely on the values for the parameters estimated 
by the mathematical technique developed in this paper. 

Addition of an apoptotic agent or withdrawal of a vital agent are expected to induce apop- 
tosis through different pathways, respectively by a caspase-initiated cascade llOj or by 
the mitochondria/cytochrome c pathway We note, however, that we cannot compare the 
concentrations of different reagents/proteins since each one of them acts in a different way in 
the cells. Different molecular processes entailing different cellular mechanisms can nevertheless 
reflect the same populational effect, thus rendering themselves to be described by the same 
mathematical model. Though at the level of description of our first-order model the parame- 
ters estimated for each mechanism are not comparable, they are nonetheless valid to suggest a 
contribution of OS-9 to cell viability and resistance to apoptosis. 

The process of explaining a mechanism involves an active process of simplification: it is nec- 
essary to make assumptions to define the most relevant components contributing to the overall 
observed dynamics, while removing factors which contribute less. The isolation and identifica- 
tion of the core principles governing the dynamics is paramount to the detailed understanding 
of any system. [For a discussion, see Ref. Following this approach, we present a minimal 

model which is capable of capturing the structure of the experimental data. This model does 
not fit perfectly the experimental data but it is built upon biologically plausible and consistent 
assumptions. The values for all the parameters in the model are calculated using the least mean 
square method and are, in that sense, optimal. More complex models could be derived which, 
by using more parameters, would provide a better fit to the data. However, such models would 
be less informative since the functional role of the additional parameters would be speculative. 
Moreover, experimental data are always noisy and a perfect fit to a specific data set means loss 
in generality. It is important to emphasise that very little is known about the dynamics asso- 
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dated with OS-9. Our minimal model conciliates the scarce information available with robust 
assumptions to produce a simple but general model capable of quantifying the contribution of 
OS-9 as a potent anti-apoptotic factor. This is a first-order model which serves as a valuable 
starting point for experiments aiming at a more detailed understanding of the role of OS-9. To 
the best of our knowledge, this is the only available model for the quantification of the OS-9 
contribution to apoptosis. 

This model can be applied to other cellular systems since it allows for the discrimination 
of effects at different stages in a complex chemical mechanism and can easily be upgraded in 
accordance to the functional information available. 
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V. CAPTIONS 



Figures denoted by (a) reproduce the data reported in Ref. p], the points depicted giving 
the readings for the different cell populations as follows: HA-expressed cell population (black), 
17V cell line (dark gray), 17VIII cell line (light gray). Figures denoted by (b) show the 
solutions for either Cy{t) or Ca{t). 

Fig. [It Number of viable cells after withdrawal of IL-3, in the presence and in 
the absence of OS-9, from FACS analysis. In (a) we reproduce the data from Fig. 2B in 
Ref. [gI. In (b) we show the solutions for Cy{t) with a = 0.057. The black curve depicts Cj,{t) 
when P = 0; the gray curves depict Cy(t), with the dark gray curve for Pnv = 0.17 and the light 
grey curve for Pnviii = 0.31. (The upper dashed curves represent the fittings to the maxima 
data with a = 0.056, Puv = 0.19 and Pnviii = 0.32; the lower dashed curves represent the 
fittings to the minima data with ao = 0.058, Pnv = 0.15 and Pnviii = 0.30.) 

Fig. [2] Number of apoptotic cells after withdrawal of IL-3, in the presence and in 
the absence of OS-9, from FACS analysis. In (a) we reproduce the data from Fig. 2C 
in Ref. Q]. In (b) we show the solutions for Ca{t) with a = 0.057 and 7 = 0.30. The black 
curve depicts Ca(t) when = 0; the gray curves depict Ca(t), with the dark gray curve for 
Pnv = 0.17 and 6nv = 0.099, and the light gray curve for Pnviii = 0.31 and Snviii = 0.086. 
(The upper dashed curves represent the fittings to the maxima data with a = 0.058, 7 = 0.26, 
l3i7v = 0.15 and 6nv = 0.070, and Pnviii = 0.30 and Snviii = 0.029; the lower dashed curves 
represent the fittings to the minima data with a = 0.056, 7 = 0.34, jSnv = 0.19 and Suv = 0.14, 
and Pnviii = 0.32 and Suviii = 0.16.) 

Fig. [3] Number of viable cells after treatment with staurosporine l/iM, in the pres- 
ence and in the absence of OS9, from FACS analysis. In (a) we reproduce the data from 
Fig. 5B in Ref. 6|. In (b) we show the solutions for c^,(t) with a = 0.20. The black curve depicts 
Cy(t) when (3 = 0. The gray curves depict Cy(t), with the dark gray curve for Puy = 0.18 and 
the light gray curve for Pnviii = 0.24. (The upper dashed curves represent the fittings to the 
maxima data with a = 0.19, Pnv = 0.14 and Pnviii = 0.23; the lower dashed curves represent 
the fittings to the minima data with a = 0.20, Pnv = 0.17 and Pnviii = 0.20.) 

Fig. [4] Number of apoptotic cells after treatment with staurosporine 1/iM, in the 
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presence and in the absence of OS-9, from FACS analysis. In (a) we reproduce the data 
from Fig. 5C in Ref. |6|]. In (b) we show the solutions for Ca{t) with a = 0.20 and 7 = 0.72. The 
black curve depicts Ca{t) when P,6 = 0; the gray curves depict Ca{t), with the dark gray curve 
for Pny = 0.18 and = 0.19, and the light gray curve for Pijviii = 0.24 and Suviii = 0.27. 
(The upper dashed curves represent the fittings to the maxima data with a = 0.20, 7 = 0.62, 
Pi7V = 0.17 and 6nv = 0.09, and Pnviii = 0.20 and Snviii = 0.20; the lower dashed curves 
represent the fittings to the minima data with a = 0.19, 7 = 0.90, Pnv = 0.14 and 6ijv = 0.33, 
and Pnviii = 0.23 and Snviii = 0.39.) 

Fig. [5] Number of viable cells after treatment with staurosporine O.l/xM, in the 
presence and in the absence of OS9, from FACS analysis. In (a) we reproduce the 
data from Fig. 5B in Ref. IgJ. In (b) we show the solutions for Ct,(t) with a = 0.15. The black 
curve depicts Cv(t) when P = 0; the gray curves depict c„(t), with the dark gray curve giving 
Pnv = 0.19 and the light gray curve Pnviii = 0.35. (The upper dashed curves represent the 
fittings to the maxima data with a = 0.15, (3nv = 0.18 and Pnviii = 0.35; the lower dashed 
curves represent the fittings to the minima data with a = 0.16, Pnv = 0.20 and Pnviii = 0.34.) 

Fig. [6] Number of apoptotic cells after treatment with staurosporine 0.1/iM, in the 
presence and in the absence of OS-9, from FACS analysis. In (a) we reproduce the 
data from Fig. 5C in Ref. 6|]. In (b) we show the solutions for Ca{t) with a = 0.15 and 7 = 0.20. 
The black curve depicts Ca(t) when (3,6 = 0; the gray curves depict Ca(t), with Pnv = 0.19 
and 6nv = for the dark gray curve, with Pnviii = 0.35 and Siyviii = for the light gray 
curve. (The upper dashed curves represent the fittings to the maxima data with a = 0.16, 
7 = 0.19, Pnv = 0.20 and 6nv = 0, and Pirviii = 0.34 and Snviii = 0; the lower dashed 
curves represent the fittings to the minima data with a = 0.15, 7 = 0.20, Pnv = 0.18 and 
Si7v = 0, and Pnviii = 0.35 and dnviii = 0.) 

Fig. [7] Number of viable cells after treatment with TNFa, in the presence and in 
the absence of OS-9, from the FACS analysis. In (a) we reproduce the data from Fig. 6B 
in Ref. 6|]. In (b) we show the solutions for c^(t) with a = 0.32. The black curve depicts Cy{t) 
when /3 = 0; the gray curves depict Ct,(t), with the dark gray curve for jSnv = 0.46 and the light 
grey curve for Pnviii = 0.54. (The upper dashed curves represent the fittings to the maxima 
data with ao = 0.31, Pnv = 0.47 and Pnviii = 0.54; the lower dashed curves represent the 
fittings to the minima data with ao = 0.33, Pnv = 0.46 and Pnviii = 0.54.) 
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Fig. [5] Number of apoptotic cells after treatment with TNFa, in the presence and 
in the absence of OS-9, from the FACS analysis. In (a) we reproduce the data from 
Fig. 6C in Ref. j^. In (b) we show the solutions for Ca{t) with a = 0.32 and 7 = 0.045. 
The black curve depicts Cait) when (3,6 = 0; the gray curves depict Ca(t), with the dark 
gray curve for Pny = 0.46 and 6nv = 0, and the light gray curve for Pnviii = 0.54 and 
617VI11 = 0. (The upper dashed curves represent the fittings to the maxima data with a = 0.33, 
7 = 0.036, Pnv = 0.46 and 6ijv = 0, and Pnviii = 0.54 and Snviii = 0; the lower dashed 
curves represent the fittings to the minima data with a = 0.31, 7 = 0.055, /Snv = 0.47 and 
6nv = 0.33, and Pnviii = 0.54 and Snviii = 0.) 

Table HI Kinetical model parameters estimated for both cell line populations. We 

summarize the values of the parameters estimated for both cell line populations and for the 
three different mechanisms of apoptosis induction. 
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VI. FIGURES 





FIG. 2: Number of apoptotic cells after withdrawal of IL-3, in the presence and in the 
absence of OS- 9, from FACS analysis. 
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FIG. 3: Number of viable cells after treatment with staurosporine 1/xM, in the presence 
and in the absence of OS9, from FACS analysis. 




FIG. 4: Number of apoptotic cells after treatment with staurosporine 1/xM, in the presence 
and in the absence of OS-9, from FACS analysis. 



22 



(a) 




5 10 15 20 25 5 10 15 20 25 



t [hours] t [hours] 



FIG. 5: Number of viable cells after treatment with staurosporine 0.1/xM, in the presence 
and in the absence of OS9, from FACS analysis. 




FIG. 6: Number of apoptotic cells after treatment with staurosporine 0.1/xM, in the 
presence and in the absence of OS-9, from FACS analysis. 
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FIG. 7: Number of viable cells after treatment with TNFa, in the presence and in the 
absence of OS-9, from the FACS analysis. 
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FIG. 8: Number of apoptotic cells after treatment with TNFa, in the presence and in 
the absence of OS-9, from the FACS analysis. 
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